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Abstract Shannon gave a lower bound in 1959 on the binary rate of spherical 
codes of given minimum Euclidean distance p. Using nonconstructive codes 
over a finite alphabet, we give a lower bound that is weaker but very close 
for small values of p. The construction is based on the Yaglom map combined 
with some finite sphere packings obtained from nonconstructive codes for the 
Euclidean metric. Concatenating geometric codes meeting the TVZ bound 
with a Lee metric BCH code over GF{p), we obtain spherical codes that are 
polynomial time constructible. Their parameters outperform those obtained 
by Lachaud and Stern in 1994. At very high rate they are above 98 per cent 
of the Shannon bound. 

Keywords spherical codes • codes for the Euclidean metric • saddle point 
method 



1 Introduction 

A spherical code is a finite set of points of the unit sphere in a Euclidean 
space of finite dimension. For motivation and background see l5l[2]. Let X 
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denote a spherical code of M". Denote by p its Euclidean squared minimum 
distance. Define its binary rate as 

R^p) :^lin,supi^^i«. 

n 

Chabauty in 1953 and Shannon in 1959 0, gave a lower bound 

R{p) > Rsip) 1 - (1/2) log2(p(4 - p)). 

Lachaud and Stern, in 1994 gave a lower bound on the rate R*{p) of poly- 
nomial time constructible spherical codes as 

R,{p) > 0.5Rs{p). 

In the present, work we shall give a lower bound based on nonconstructive 
lattice packings 

R{p)>Rl{p) :=-(l/2)log2(p). 
It can be shown by direct substitution that 

Rsip) > Rl{p) 

for ah < p < 4. 

However, for small p, the two curves are very close to each other; 

R{p)^Rl{p)=0{p^). 

We give a family of spherical codes based on non constructive codes for 
the Euclidean metric whose asymptotic performance is just as good. 

R{p) > RUp), 

This shows that using finite alphabet codes for constructing finite packings 
is as efficient as using truncation of dense infinite sphere packings. 

However, the main result of this work is a lower bound based on an explicit 
family of polynomial time constructible codes for the Euclidean metric that 
outperforms the Lachaud Stern bound almost by a factor of 2 

R,{p) > 0.98i?s(p). 

for some very small values of p in the range < p < 1. These codes, 
furthermore admit efficient encoding and decoding algorithms. 

The article is organized as follows. In the next section we explore the 
Yaglom map and use it on truncations of infinite lattices. In Section HI we 
study the performance of long codes for the Euclidean metric over a fixed 
alphabet. The analysis technique is based on the saddle point method. In 
Section IV we give the explicit construction of long such codes and prove it is 
better not only than the relevant Gilbert bound but also after Yaglom map 
than 98 per cent of the Shannon bound. 
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2 The Yaglom map 

Following Yaglom [2J Chap. 9, Thm. 6], we inject M" into M"+^ by the map 

Y : X 1-^ {x, ^/R^ — X ■ x). 
Note that this map sends the ball B{n, R) of radius R that is 

B{n,R) := {x G R"\x- x < i?^} 

in K" into the sphere S{n, R) of radius R in M"+-'^ 

S{n,R) ■.= {xe«"+^\x-x^R^} 

and that distance between points can only increase. The following result is 
then immediate. 

Proposition 1 If P is a packing of spheres of diameter d in B{n,R) then 
Y{P) is a spherical code in S(n,R) of minimum square distance d^. 
Using lattice packings we can prove our first bound. 

Theorem 1 There are families of spherical codes built from lattice packings 
by the Yaglom map such that 

R{p)>Rl{p) :=-(l/2)log2(p). 



Proof. By the Yaglom bound Chap. 9, Thm. 6] we know there are spherical 
codes satisfying 

R>6 + 1 + Rl{p), 

where 6 is the asymptotic exponent of the best density of lattice packings in 
E", or in the notation of jSj, Chap. 1]. Now by the Minkowski bound 

Chap. 1, (46)], we know that S> -1. The result follows. □ 
Using lattice packings built from codes over fields by the so-called Con- 
struction A [21 p. 182] we obtain our second bound. 

Theorem 2 There are families of spherical codes built from Construction A 
lattice packings such that 

R{p)>Rl{p). 

Proof. 

By [12j we know there are lattice packings built from codes over prime 
fields by Construction A having asymptotic density 

5 > -1. 

The result follows, as the preceding one, upon applying the Yaglom bound. 

□ 
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This result will be obtained also for finite sphere packings made of codes 
for the Euclidean metric (Theorem ([SJ). Using explicit lattice packings con- 
structed by AG techniques we can prove our last bound. Note that, even 
though these lattices are polynomial time constructible, the resulting spheri- 
cal codes may not be because of the cost of truncation which may require an 
exponential time algorithm like Finke-Pohst [5]. 

Theorem 3 There are families of spherical codes of binary rate R{p) built 
from lattice packings by the Yaglom map such that 

>i?L(p)-1.30. 

Proof. By the Yaglom bound [U Chap. 9, Thm. 6] we know there are spherical 
codes satisfying 

R>S + 1 + Rl{p), 

where S > —2.30 is the asymptotic exponent of the best density of lattice 
packings constructed in [S]. The result follows. □ 

Corollary 1 For p < e~^'^ there are spherical codes constructed from polyno- 
mial time constructible lattices with R > 0.98Rs{p)- 

Proof. We apply the preceding Theorem with the approximation Rs{p) ~ 
Rl{p)- Wc obtain the said bound for p < 2"^^° w e~^°. □ 

3 Codes for the Euclidean metric 

We consider codes of length n over the integers modulo an integer q. If q = 
2s -I- 1 is odd we represent Zg on the real line by the constellation 

{-s,--- ,-1,0,1, ••• ,s.} 

If q = 2s + 2 is even we represent on the real line by the constellation 

{-s - 1/2, • • • , -3/2, -1/2, 1/2, • • • , s - 1/2, s + 1/2}, 

which is a shift of the natural representation 

{-s,--- ,-l,0,l,--- ,s,s-M} 

by —1/2. We denote by (j> the induced map from Z^ into M". If q is odd then 
(l){x) ■ (j){x) < ns^, while if q is even <j)[x) ■ (j){x) < n{s + 1/2)^, We consider the 
Euclidean distance dsi,) on Z^ defined by dE{x,y) = i't'i^) ~ 4>{y)Y ■ This 
distance could have been defined as induced by the standard Euclidean weight 
we{) on Zg that is we{x) — min(a;^, {q — x)^. 

Proposition 2 If C is a code of length n and minimum Euclidean distance 
d over Iq then (j){C) is a packing of spheres of diameter y/d of Bin, Sy/n) if q 
is odd and of B{n, (s -I- 1/2)y^) if q is even. 
Combining with Proposition [1] we obtain 
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Proposition 3 If C is a code of length n and minimum Euclidean distance 
d over Iq then Y{<f>(C)) is a spherical code of squared Euclidean distance d of 
S{n,Sy/n) if q is odd and of S(ri, {s + l/2)y/n) if q is even. 

Let V{n, q, r) denote the size of the bah of radius r for the Euclidean metric 
in Zq. The following result is the analogue of the standard Gilbert bound in 
that setting. 

Proposition 4 There are codes in Tll^ of Euclidean distance d and of cardi- 
nality 



\C\> 



V{n,q,d-1)- 



The technically difficult part is to estimate the asymptotic exponent of 
V{n,q,r) for n large and r = [AnJ. To that end, we introduce the following 
generating series f{z) = 1 + 2s{z) with 



for odd g and 



z 

i=l 



s-l 

Z*' + z' 
i=l 



for even q. 

Theorem 4 The asymptotic exponent of V{n, q,r) for n large and r = [AnJ 



IS 



,. Aog2{V {n,q,r)) /./ \, / n 

lim( ) = log2(/(/i)) - Alog2(Ai) 

where fi is the unique real positive solution of zf'{z) = Xf(z). 
Proof. The generating function for the numbers V(n,q,r) is of the form 
f{z)"g{z) with g{z) = 1/(1 — z). The result follows by application of [U Corol- 
lary 2]. □ 

Theorem 5 With the preceding notation, there are spherical codes of the unit 
sphere S{n + 1, 1) with relative squared Euclidean distance p < 1 of binary rate 

R > log2(q) - log2(/(Ai)) + aplog2(M) 

where (i is the unique real positive solution of zf'(z) ~ apf{z), with a = 
for q odd and a = {s + 1/2)^ for q even. 

Proof. We combine the Gilbert bound of Proposition U with the estimate 
of Theorem [4] to construct long Euclidean metric codes with prescribed pa- 
rameters and the Yaglom map of Proposition [3] to derive spherical codes from 
them. Note that we let A = ap to rescale S{n + 1, ^/na) to S{n + 1, 1). □ 
This approach can be generalized to the case of a large varying alphabet. 
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Theorem 6 There are spherical codes, constructed from codes over Z^, with 
q odd and variable and asymptotic rate 

R{p) > Rl{p)~c, 

with cKi 0.77 X 10"^. 

Proof. We use the same saddle point analysis as [TT] on the generating 
function 

QO 

1 + 2^/. 

1=1 

Note that s does not occur in the latter expression since q > 2r + 1. □ 



4 Constructive bound 

By the TVZ bound [3 Th. 13.5.4] we know there are families of geometric codes 
over GF{Q) for Q a square with rate TZ and relative distance A satisfying 

7^ + zi > 1 - ^ 



- 1 

We concatenate this geometric code with a code over Zq of parameters 
[n, fc] and minimum Euclidean distance dE- We must assume therefore that 
Q = q*^. If g is not a prime, we label 1q by the elements of GF{q) in an 
arbitrary fashion. In order to apply the TVZ bound we must assume Q to be 
square. 

Proposition 5 With the above notation, the Yaglom map of the concatenated 
code has parameters {R, p) above the straight line 

Rn pns"^ 1 



\og2{q)k dE 

Proof. We assume some familiarity with concatenation [71 §5.5]. The g— ary 
rate of the concatenated code is ^TZ. The binary rate of the spherical code is 
therefore R = log2(g)^72,. The relative distance of the inner Euclidean metric 
code is The relative distance of the concatenated code is After 
normalization to reduce to the unit sphere Proposition ([2]) yields 

P > 

Substituting into the TVZ bound we are done. □ 
Let us take, as inner codes, the Roth-Siegel BCH codes [TU]. The Lee 
minimum distance is an immediate lower bound on the Euclidean minimum 
distance of such a code. 
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Theorem 7 For each prime p > 7 and every integer 1 < t < {p + l)/2, such 
that p is congruent to t + 1 mod 2, there is a family of spherical codes with 
binary rate R and minimum squared Euclidean distance p over Zj, satisfying 

R{p-l) , p(p-lf ^ ^ 1 



(p-t- l)log2(p) 8i - p(p-t-l)/2 _ i' 

Proof. For these codes n ^ p — 1, >2t and fc > n — i, by [TU]. □ 
We need the foUowing elementary result whose proof is omitted. 

Proposition 6 The equation of the tangent in p of the curve {p, XRl{p)) is 
of the form 

X Y ^ 

with 

A^p{i- Hp)) 

We are now in a position to state and prove the main result of this note. 

Theorem 8 Take p = 5432455719452623343140299649993224712642268405087 
972148236533041723675544652674874508958455203602044198462638584629866 
4106668659730094751 and t/n = : 00155359 in Theorem For all values 
of P l£ e~640 : 48 and A ~ 0.98 the bound of Theorem ^ is strictly above the 
tangent to the curve (p, XRl{p)) at p. 

Proof. Denote by / = 1 — pip-t-i)/2_i ^ number very close to 1. Put t ^ t/n. 
We find values of p and p that the intersections of the tangent with the axes 
are less than 

(p-i-l)log2(p) 



A < 



p~l 
Srf 



Approximating p — 1 hy p and / by 1 we get by replacing A, B by their 
expressions from Proposition |6l 

A(l-ln(p)) <2(l-r)ln(p) 
8t 

p{i - Hp)) < — 

Write 

X = ln(p) 
y = ln(p) 

After this change of variables, and getting rid of r, between 
(1 - x) exp(a; + 2y)/8 < r < 1 - A(l - x) /2y 
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wc obtain the equation 

exp(a; + 2y)(l - x) + 4A(1 -x)/y< 8, 

which can be solved numerically and graphically as shown in Figure [T] Once 
the values of x and y have been constrained in that way the assertions of the 
Theorem are checked by direct (machine) computation. □ 

315.10 P ' ' ' ' ^ ' ' ' ' ^ ' ' ' ' ^ ' ' ' ' ^ 



315.05 




314.80 h I , , , , I , , , , I - 

-640.8 -640.7 -640.6 -640.5 -640.4 

ln(p) 

Fig. 1 Attainable Region for Inp and Inp; here, A = 0.98 

Recall that the envelope of a family of curves is a curve all the points of 
which are tangent at some point to one of the curves in the family. We can 
obtain the envelope of the family of straight lines promised by TVZ when p 
varies as follows 

Theorem 9 Some points produced by Theorem lie on the curve 

(l-T)y l-X 

where 
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X = ln(p) 

y = ln(p) 

X + 2y = c 

8t = (1 — x) exp(c) 



and c is a constant depending on the range of t. Proof. Keep the notation of 
the preceding proof. An approximate equation for the straight hnc of Theorem 
His 



(l-r)ln2(p 8r 
We choose arbitrarily 

X + 2y ~ c 

for some constant c and 

8r = (1 — x) exp(c). 
The result follows afer some algebra. 




Experimentally, the curve obtained in Figure [5] shows that our constructive 
codes outperform Ax the Shannon bound. Note that we took, here, a value of 
A = 0.976 to show more clearly the performance of such codes, even if these 
codes have been designed for A = 0.98. 
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5 Complexity issues 

While the size of the large prime in Theorem 8 might seem outlandish, we 
believe that our constructions of spherical codes might have some practical 
interest for more realistic inner codes. We observe that the complexity of the 
Yaglom map is 0{N) and that the concatenated codes being linear, encoding 
of the spherical codes is linear in time and space. This contrasts with most 
standard constructions of spherical codes, especially those based on constant 
weight codes [S]. As for decoding, the decoding of concatenated code is possible 
knowing decoding algorithms for the inner and outer codes [3| . As explained in 
[7j §13.5.2] polynomial time decoding algorithms are known for certain codes 
obtained from the Garcia Stichtenoth towers of function fields. An algebraic 
decoding algorithm up to the Lee error correcting capacity is provided for the 
BCH codes we use in [TU] . 



6 Conclusion 

In this work we have constructed spherical codes by using finite packings of 
spheres within a ball and mapping them on the surface of the sphere in the 
next dimension by using the Yaglom map. The construction of finite packing 
was done by using codes for the Euclidean metric. We hope this will stimulate 
research in this area. For instance can the minimum Lee distance bounds of 
[TU] for BCH codes be improved for the Euclidean metric? Are there perfect 
codes for the Euclidean metric? This might provide stronger inner codes for the 
concatenation construction and might lead to an improvement of the Shannon 
lower bound. 

Acknowledgement The authors thank Christine Bachoc and Philippe 
Gaborit for helpful discussions. 
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